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Abstract 

We investigate oscillating instanton solutions of a self-gravitating scalar 
field between degenerate vacua. We show that there exist 0(4)- 
symmetric oscillating solutions in a de Sitter background. The geometry 
of this solution is finite and preserves the Z2 symmetry. The nontrivial 
solution corresponding to tunneling is possible only if the effect of gravity 
is taken into account. We present numerical solutions of this instanton, 
including the phase diagram of solutions in terms of the parameters of 
the present work and the variation of energy densities. Our solutions 
can be interpreted as solutions describing an instanton-induced domain 
wall or branewor Id-like object rather than a kink- induced domain wall 
or braneworld. The oscillating instanton solutions have a thick wall and 
the solutions can be interpreted as a mechanism providing nucleation 
of the thick wall for topological inflation. We remark that Z2 invariant 
solutions also exist in a flat and anti-de Sitter background, though the 
physical significance is not clear. 
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1 Introduction 



The interpolating solution between degenerate vacua in curved space can be the solution repre- 
senting formation of a domain wall or a braneworld-like object with Z2 symmetry if the thin- wall 
approximation scheme is used in the theory. Moreover, the solution after the analytic continu- 
ation expands, from an observer's point of view on the wall, without eating up bulk (inside and 
outside) spacetime. Can the braneworld or domain wall have a nontrivial internal structure? 
The structure can be made with oscillating instanton solutions. What is the meaning of this 
internal structure? Can the number of oscillations indicate certain criteria based on specific 
parameters? In order to get the answers to the above questions, we need to obtain oscillating 
solutions, interpret the meaning of the solutions, and construct a phase space of all allowed 
oscillating solutions. 

In the absence of gravity, instantons are usually defined as solutions with a finite action 
to the classical field equations in Euclidean space obeying appropriate boundary conditions. 
The Yang-Mills instantons are also finite action solutions and have been studied extensively in 
gauge theories [Tl|2] as well as in string theory [3lll], and references therein. In gravitational 
theory, there are several kinds of instantons [5]. Instantons should also be solutions to the 
Euclidean Einstein equations. Some instantons usually cannot be analytically continued to 
Lorentzian. Others can be analytically continued to Lorentzian. Instantons may be related to 
the semiclassical description of quantum gravity. However, this issue is beyond the scope of the 
present work and is therefore not being discussed here. 

In this work, we consider tunneling phenomena in a double-well potential. Classically, a 
particle trapped in one well cannot penetrate through the potential barrier of that well, thus 
unaffected by the presence of the other well. However, if the central potential barrier is not 
infinitely high, there is tunneling between the two minima. The tunneling process is quantum 
mechanically described by the Euclidean solution obeying appropriate boundary conditions. 
The Euclidean solution interpolates between two different classical vacua. There exist two 
kinds of Euclidean solutions describing tunneling phenomena in the double-well potential. One, 
which is for tunneling in an asymmetric double-well potential, corresponds to a bounce solution. 
The bounce solution describes the decay of a background vacuum state. The other, which 
is for tunneling in a symmetric double-well potential, corresponds to an instanton solution. 
The instanton solution, in this case, describes a general shift in the ground state energy of 
the classical vacuum due to the presence of an additional potential well, lifting the classical 
degeneracy [6]. 

The instanton solution in one dimension is equivalent to a static soliton in (1 + 1) dimen- 
sions. Thus, the method for kink solution in two dimensions can be employed for studying the 
instanton solution in one dimension. The Euclidean equation of motion can be treated like a 
one-particle equation of motion with an evolution parameter playing the role of time in the 
inverted potential. A particle starts at the top of one hill at minus infinity in Euclidean time 
and arrives at the top of the other hill at plus infinity. The Euclidean action for the tunneling 
solution can be calculated easily using the fact that the Euclidean energy of an instanton is 
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equal to zero in the system. This action is identical to the energy of a static solution of the 
(1 + l)-dimensional soliton theory. In addition, the action has the same value as the action ob- 
tained in connection with the WKB approximation of the splitting in energies of the two lowest 
levels for the symmetric double-well potential. In semiclassical approximation, the action is 
dominated by classical configuration in evaluating the path integral. 

For field theoretical solutions in (3 + 1) dimensions, the Euclidean equation of motion for 
0(4) symmetry has an additional term, which can be interpreted as a damping term. Thus, 
it is difficult to obtain the 0(4)-symmetric instanton solution in this system when gravity is 
switched off. The damping term can be changed into the antidamping term not only in de 
Sitter (dS) space but also in both flat and anti-de Sitter (AdS) space when gravity is taken 
into account [7]. The solutions have exact Zi symmetry and give rise to geometry of a finite 
size. Using a homogeneous and static scalar field, which has a constant value of $ everywhere 
in Euclidean space, was shown to be impossible for the solution to have an infinite size due to 
the infinite cost of energy. 

The bounce solution is related to the nucleation of a true (false) vacuum bubble describing 
decay of a background vacuum state. The process has been studied within various contexts for 
several decades. It was first investigated in Ref. [8] and developed in both fiat [9] and curved 
spacetime [TUIIII]. A homogeneous Euclidean configuration in which the scalar field jumps 
simultaneously onto the top of the potential barrier was investigated in Ref. [12] and further 
studied in Ref. [13]. As a special case of the true vacuum bubble, a vacuum bubble with a finite- 
sized background after nucleation was studied in Ref. [H]. The decay of false monopoles with 
a gauge group was also studied using the thin- wall approximation [15]. The bubble or brane 
resulting from fiux tunneling was studied in a six- dimensional Einstein-Maxwell theory |16] . 

The mechanism for nucleation of a false vacuum bubble in a true vacuum background has 
also been studied within various contexts. Nucleation of a large false vacuum bubble in dS space 
was obtained in Ref. [TT] and nucleation with a global monopole in Ref. [TS]. The mechanism 
for nucleation of a small false vacuum bubble was obtained in the Einstein gravity with a 
nonminimally coupled scalar field [19], with Gauss-Bonnet term in Ref. [20], and using Brans- 
Dicke type theory [21]. The classification of vacuum bubbles including false vacuum bubbles 
in the dS background in the Einstein gravity was obtained in Ref. [22], in which the transition 
rate and the size of the instanton solution were evaluated in the space, as the limiting case of 
large true vacuum bubble or large false vacuum bubble. 

The oscillating solution with 0(4) symmetry in dS space was first studied in Ref. [23], where 
the authors found the solution to oscillating scalar field $ in fixed background. They adopted 
the fact that the role of damping term in the particle analogy picture can be changed into that 
of antidamping term in dS space if the evolution parameter exceeds half of a given range. The 
oscillation means that the field in their solutions oscillates back and forth between the two sides 
of the potential barrier. To obtain the general solutions including the effect of the backreaction, 
we will solve the coupled equations for the gravity and the scalar field simultaneously. 

In the absence of gravity, the spectrum of small perturbations about the bounce solution 
has exactly one negative mode [23]. When gravity is taken into account, the problem becomes 
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more complex [25]. On the other hand, the spectrum about the instanton solution between 
degenerate vacua has no negative mode. Thus we expect that the spectrum about our solutions 
will not have any negative modes. 

The paper is organized as follows: in the next section we set up the basic framework for this 
work. We consider the tunneling process for the symmetric double-well potential and study 
the boundary conditions for our present work in detail. The boundary conditions analyzed in 
this work can be applicable to other cases including bounce solutions. In Sec. Ill, we present 
numerical oscillating instanton solutions by solving the coupled equations for the metric and the 
scalar field simultaneously. We show that there exist 0(4)-symmetric oscillating solutions in dS 
background. We also examine the variation in the thickness of the wall and the energy density 
as the number of oscillations increased. Additionally, we examine numerically which solutions 
among oscillating solutions are relatively probable. The numerical solution is also possible in 
fiat and AdS space as long as the local maximum value of the potential is positive. In Sec. 
IV, we analyze the properties of oscillating instanton solutions and construct the phase space 
of all our solutions in terms of the two parameters for this work. The phase space of solutions 
exhibits the behaviors that occur when the number of oscillations increases and indicates the 
regions where there are no solutions depending on the parameters. In Sec. V, we summarize 
and discuss our results. We interpret the physical meaning of our oscillating solutions with 
thick walls and discuss the probability of our solutions. 



2 The setup and the boundary conditions 

The vacuum-to-vacuum transition amplitude called the generating functional for the Green's 
function can serve as a starting point for the nonperturbative treatment of the theory. In the 
path integral formalism, it can be visualized as the summation over all possible paths moving 
from the initial to the final state. The transition amplitude in the semiclassical approximation 
is dominated by those paths for which the Euclidean action difference is stationary. Thus, the 
tunneling amplitude can be evaluated in terms of the classical configuration and represented 
as Ae~^^ in this approximation, where the leading semiclassical exponent AS is the difference 
between the Euclidean action corresponding to a classical solution and the background action 
itself. The prefactor A comes from the first order quantum correction [ni[21] • 
Let us consider the following action: 



J M 



Zk z 



V-hd'^x 

dM 



3 K-Ko 



where g = detg^^, k = SnG, R denotes the scalar curvature of the spacetime Ai, K and Ko are 
the traces of the extrinsic curvatures of dM. in the metric g^j^i, and 77^,^, respectively, and the 
second term on the right-hand side is the boundary term [26]. The gravitational field equations 
can be obtained properly from a variational principle with this boundary term. This term is 
also necessary to obtain the correct action. 
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We reconsider the tunneling process for the symmetric double-well potential in curved space 
similar to our previous work [7]. The potential U{^), which represents the energy density of a 
homogeneous and static scalar field, has two degenerate minima 



A 



A 



Uo. 



(2) 



The cosmological constant is given by A = kUo, hence the space will be dS, fiat, or AdS 
depending on whether Uq > 0,Uo = 0, or Uq < 0. We will make an attempt to obtain oscillating 
instanton solutions describing tunneling between degenerate vacua in these backgrounds. 

To evaluate AS* and show the existence of the solution, one has to take the analytic con- 
tinuation to Euclidean space. We assume the 0(4) symmetry for both the geometry and the 
scalar field as in Ref. (TU] 



ds^ = drf + p2(ry) [dx^ + sin^ x [dO^ + sin^ Odcl)'^)] . 



(3) 



In this case $ and p depend only on rj, and the Euclidean field equations for them can be 
written in the form: 



$' + — $' = — and p" = - 
p d^ 

respectively and the Hamiltonian constraint is given by 

n2 



K 



'-p{^'' + U), 



p 



12 



1 - 



Kp 



2 



U 



0. 



(4) 



(5) 



In order to yield the solution the constraint requires a delicate balance among the terms. 
Otherwise the solution can provide qualitatively incorrect behavior [27] . 

Now we have to consider the boundary conditions to solve Eqs. ([4]) and ([5]). There are two 
different methods. The first one is an initial value problem, in which we impose the initial 
conditions for the equations. For this to work, initial conditions are provided for the values of 
the fields p and $ or their derivatives p' and $' at ?7 = as follows: 



p\r,=o = 0, 



dp 
df] 



1, $ 



ri=0 



r)=0 



$n, and 



dr] 



0, 



ri=0 



(6) 



where the first condition means that the space including a solution is a geodesically complete 
space. The second condition stems from Eq. ([5]). The fourth condition $' = is due to 
the regularity condition as can be seen from the first equation in Eq. (j4]). However, the third 
condition for the initial value of $ is not determined. One should find the initial value of $ using 
the undershoot-overshoot procedure. This procedure particularly useful for the nonoscillating 
bounce solution, i.e. a one-crossing solution as we will explain about it in the next section. 
Some initial $0 will give the overshooting, in which the value of $ at late 77 value will go beyond 
fi/^/X. Some other initial value of $0 will give the undershoot, in which the value of $ at late 
r] does not reach jj/y/X. Thus, the value $0 must be within an intermediate position between 
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the undershoot and overshoot [9]. The existence of oscillating solutions is hidden within the 
undershoot as we will see in the next sections. The other method is a boundary value problem 
or a two boundary value problem. We impose conditions specified at 77 = and rj = rj^ax- For 
this to work, we choose the values of the field p and derivatives of the field $ as follows: 



p\n=o — 0, p\n=rimaa: " — 



dr] 



d^ 

= 0, and — 

ri=o dr] 



= 0. (7) 



'Umax is the maximum value of 77 and will have a finite value. These conditions are useful for 
obtaining solutions with Z2 symmetry. 

In order to solve the Euclidean field Eqs. (jl]) and ^ numerically, we first rewrite the 
equations in terms of dimensionless variables as in Ref [19] : 



^ = ;7($), -^ = ^\ PV = V, I^P = P, and = k. (8) 

p p, A 

These variables give 

f/(<l) = ^f<i>2-l)' + [/o, (9) 



and the Euclidean field equations for $ and p are rewritten as: 

3p'~, _ dU K^,~,2 



$" + = ^ and p" = p($'' + U), (10) 

p d^ 3 



respectively. The Hamiltonian constraint is given by 

p--l-^Q$--t/)=0. (11) 

We now make comments on the boundary conditions more precisely to solve Eqs. (fTOl) and 
( JTTl) numerically. As we already mentioned in the boundary value problem, there are two kinds 
of conditions. The relaxation method is one of the methods employed to solve a two boundary 
value problem. However we do not know the exact value of fjmax in our case, thus we can not 
impose the boundary condition at 77 = fjmax- In this work, we employ the shooting method 
using the adaptive step size Runge-Kutta as an initial value problem as done in Ref. [28]. For 
this procedure we choose the initial values of ^{fjiniuai), ^' iViniUai) , piViniUai), and p' {fjiniuai) at 
V = fjinitiai as follows: 

^{f|^n^Ual) ~ l-Q + ^^^O ($0 " l) + (3$o 



$'(to/) ~ -$o($^-lj+-(3$^-l), (12) 

PiViniUal) ~ e, 
p' {fjinitiai) ~ 1, 
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where fjiniuai = + e and e ^ 1 . If we find tfie initial value $o > the other conditions are given 
by Eqs. fll2p . Furthermore we impose additional conditions implicitly. To avoid a singular 
solution at f] = fjmax in Eq. (fTUj) and to demand the Z2 symmetry, the conditions d^/dfj — 
and p — J- as — » fjmax are needed. In this work, we require that the value of d^/dfj goes to a 
value smaller than 10~^ as 77 — fjmax, as the exact value of fjmax is not known. 

We will examine the energy density of the oscillating region of the solution to observe a 
domain wall structure. The action for the one- dimensional instanton has the following form 



POD 


1$'2 + U 


J —00 


2 



dr 



2/i3 
'3X'' 



where we take the potential similar to the one in Eq. To obtain a solution with finite 
action, it is required that $' — )■ and ?7 — )■ as r — )■ ±00. This action is equal to the surface 
tension or the surface energy density of the wall [9]. When gravity is taken into account, the 
situation is more complicated. To simplify things, we only consider the Euclidean action of the 
bulk part in Eq. ([1]) to get. 



gsd'^XE 



M 



2k 2 



U 



27r2 / p'dri[-Ul 



(13) 



where Re = 6[l/p^ — p'^/p^ — p" / p] and we used Eqs. (jl]) and ([5]) to arrive at this. Thus the 
contributions coming from the geometry and kinetic energy in the Euclidean action are included 
to be the potential effectively. On the other hand, the action in Ref. [10] can be rewritten as 



127r2 



pdr] 1 



Kp^U 



(14) 



where the authors used integration by parts. Since they were only interested in the difference 
between the actions of the two solutions that agree at infinity, the surface term evolving from 
the integration by parts is regarded as unimportant. The minus sign appeared in the action is 
related to the fact that the Euclidean action for Einstein gravity is not bounded from below, 
which is known as the conformal factor problem in Euclidean quantum gravity [29]. However, 
this issue is not central to our problem and therefore we will not discuss this issue furthermore 
in the present paper. In this work, we will examine the relative probability among the solutions. 
We employ Eq. (fT3|) for the present work. The volume energy density has the following form 



Re 1 j;/9 tt 

2k 2 



U. 



(15) 



We will examine the change in the density with respect to the evolution parameter rj in the 
next section. 
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3 Oscillating instanton solutions between degenerate vacua 



In Ref. [23], the authors found an oscillating solution between the dS-dS vacuum states. They 
have shown that there exist two kinds of oscillating solutions, one is endowed with an asymmet- 
ric double-well potential while the other is with a symmetric potential. They showed how does 
the maximum allowed number Umax depend on the parameters of the theory, in which n denotes 
the crossing number of the potential barrier by the oscillating solutions. However, they solved 
the equation for the scalar field in the fixed dS background. This is a good approximation only 
for Uo ^ U{^) — Uo [23] ■ In this paper, we will solve the coupled equations for the metric 
and the scalar field simultaneously for any arbitrary vacuum energy Uo- In the first place, we 
present the details of oscillating instanton solutions including numerical solutions and the vari- 
ation of the energy density to see the structure of the domain wall and the oscillating behavior 
depending on k and kUo in dS space with Uo > 0. Then, we present oscillating solutions in flat 
space with Uo = 0. Finally, we present oscillating solutions in the AdS space with Uo < 0. 

Figure [1] shows oscillating instanton solutions representing tunneling starting from left vac- 
uum state in dS-dS degenerate vacua. The first figure illustrates the potential, in which the 
number n denotes the number of the crossing or the number of oscillations. We take Uo = 0.5 
and k = 0.04 for all the cases. The second figure illustrates the initial point $0 for each 
number of oscillations. As expected, the number of oscillations increases as the initial point, 
^iVinitiai) = ^o, moves away from the vacuum state. The third figure illustrates the solutions 

of p. The solution of p is ■sj^- sin sj'^-f] in fixed dS space. Thus, the graph of a sine type 

function near the vacuum states indicates dS space. We can see that the size of the geometry 
with such a solution decreases as the number of crossing increases because the period of the 
evolution parameter fjmax decreases as the starting point moves away from the vacuum state. 
Figure (i) illustrates the one-crossing solution, n = 1, of the field $. The one-crossing solu- 
tion corresponds to the instanton solutions with the 0(4) symmetry between the degenerate 
vacua [7], for which it was shown that the solutions exist not only in the dS space but also 
in the flat and AdS space if the local maximum value of the potential is positive. Figure (ii) 
illustrates the two-crossing solution, n = 2. The two-crossing solution, n = 2, was considered 
as a type of the double-bounce solution or anti-double-bounce solution [30], in which the au- 
thors interpreted the double-bounce solution as the spontaneous pair-creation of true vacuum 
bubbles separated by a wall. However, the solution in this paper is quite different from the 
double-bounce solution found in Ref. [30]. Since our solution of $ does not asymptotically 
approach the other vacuum state, it is difficult to interpret our solution as the double-instanton 
solution or the spontaneous pair-creation of instanton solutions. Figures (iii)-(vi) illustrate the 
n-crossing solution, for n = 3, 4, 5, 6. The maximum number n^ax = 6 depends on the param- 
eters Uo and k. In other words, there is no solution for n > rimax with Uo = 0.5 and k = 0.04. 
Figures (i), (iii), and (v) illustrate the tunneling starting from the left vacuum state to the right 
vacuum state. Figures (ii), (iv), and (vi) illustrate solutions going back to the starting point 
after oscillations. The maximum number of oscillations is determined by the parameters k and 
Uo as was observed in Ref. 1231. 
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Figure 1: The numerical solutions represent oscillating instanton solutions between dS-dS degenerate vacua. 

Figure [2] shows the variation of terms, p', p", and dU/d^, with respect to fj 

appearing in Eqs. ffTOj) and f|TT]) . In figures (i) - (vi), we see that the change of sign of p' from 
positive to negative occurs at the half period due to the Z2 symmetry. The value of p' at the 
initial and the final value of fj means there is a dS space at that point. The value of p' spans 
from —1 to 1 in all figures. The transition region of p' means the rolling duration in the inverted 
potential. In this region, all other terms have dynamical behavior. The value of $" representing 
an acceleration of the particle in the inverted potential increases, decreases, and becomes zero 
at the half period. The graph is odd function. The value of p" is always negative or zero as an 
even function according to Eq. ( ITOl) . The damping term also increases, decreases, and becomes 
zero as an odd function. The term dU/d^ has got the same property. The figures (i), (iii), 
and (v) representing tunneling show that p', damping term, and dU/d^ change their sign 
simultaneously at the half period. While the figures (ii), (iv), and (vi) representing solutions 
going back to the starting point show that only p' change its sign at the half period. All of the 
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Figure 2: Variation of terms in equations of motion between dS-dS degenerate vacua. 

behaviors represented in each figure can be well understood bearing Z2 symmetry in mind. The 
initial and final regions of p' in each figure exhibit cosine type function as the solution near the 
vacuum states indicates dS space. 

The behavior of the solutions in the ^(j])-^' (j]) plane using the phase diagram method is 
shown in Fig. [3l The figure (i) illustrates the phase diagram of a one-crossing solution, in which 
the trajectory is restricted to the upper half region of the diagram. It is the turning point 
from the damping phase to the antidamping phase when $' takes the maximum value and $ 
attains the first zero. The value of $ spans from —1 to +1 and $' from zero via maximum 
value to ~ 0.47, to zero with symmetry about the y axis. The figure (ii) illustrates the diagram 
of a two-crossing solution, in which the trajectory does not reach the opposite point $ = 1 
but return to the starting point of $ = —1. When the trajectory goes back, $' is negative 
with a symmetry about the x axis. It is the turning point from the damping phase to the 
antidamping phase when <!>' reaches the second zero and $ takes the positive value. The figure 
(iii) illustrates the diagram of a three-crossing solution. It is the turning point where $' takes 
the negative maximum value and $ attains the second zero. The figure (iv) illustrates the 
diagram of a four-crossing solution. It is the turning point when $' reaches the third zero and 
$ takes the negative value. The figure (v) illustrates the diagram of a five-crossing solution. It 
is the turning point when $' takes the positive value and $ attains the second zero. The figure 
(vi) illustrates the diagram of a six-crossing solution. It is the turning point when $' reaches 
the third zero and $ takes the positive value. The figures (i), (iii), and (v) have a symmetry 
about the y axis, whereas figures (ii), (iv), and (vi) have a symmetry about the x axis. The 
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Figure 3: The behavior of the solutions in the plane using the phase diagram method. This case belongs 
to the tunneling between dS-dS degenerate vacua. 

maximum value of $' decreases as the number of crossing increases. 

Figure S] shows the diagram of the energy density for each solutions. In each figure part, 
the solid line denotes the Euclidean energy E^^ and the dotted line denotes the volume energy 
density ^ in Eq. ( !T5|) . The Euclidean energy signifies the value after the integration of variables 
except for rj in the present case. The peaks represent a rolling phase in the valley of the inverted 
potential. The maximum value ^max is equivalent to Utop- The number of peaks is thus equal to 
the number of crossing. The peaks broaden in their range near Utop as the number of crossing 
increases. The Euclidean energy also has peaks. However, the shape of the peaks becomes 
smooth and broadens as the number of crossing increases. As can be seen from figure (vi), the 
thickness of the wall increases as the number of oscillations increases. The mountain-shaped 
graph of the Euclidean energy in each figure part is due to the integration of variables in the 
dS space. 

We now consider oscillating instanton solutions between fiat-flat degenerate vacua with 
JJo = 0. The numerical solutions in this case are shown in Fig. El where we take k = 0.2 for 
all the figures. The figure (i) illustrates the solution of $ for a one-crossing solution and the 
figure (ii) for a two-crossing solution. The final figure illustrates the solution for each p. The 
solution of p is ?7 in fixed flat space. Thus, the graph of the linear function near the vacuum 
states indicates fiat space. The maximum number rimax is 6 for the parameters Uo = and 
K = 0.2. 

Figure [6] depicts a collection of diagrams including the variation of terms in the equations of 
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Fi gure 4: The diagram for energy density of each solutions. In each figure, the dotted hne denotes the volume 
energy density ^ and the solid line denotes the Euclidean energy i?^ at constant t]. 



motion, phase diagrams, and the diagram of energy density between flat-flat degenerate vacuum 
states. Figure (i) corresponds to a one-crossing solution and figure (ii) corresponds to a two- 
crossing solution. In the second figure, the maximum value of is about 0.49. The volume 
energy density and the Euclidean energy are always positive. The initial and final region of p' 
in the first and fourth figures exhibit a horizontal nature as the solution near vacuum states 
indicates fiat space. In Figs. [5] and [6|, two solutions with n = 1 and n = 2, are illustrated and 
other solutions are omitted because the general behaviors for these cases are similar to those 
in dS space. 

Next, we consider oscillating instanton solutions between AdS-AdS degenerate vacua with 
Uo < 0. The numerical solutions in this case are shown in Fig. [TJ where we take Uo = —0.02 
and k = 0.4 for all the figures. Figure (i) illustrates the solution of $ for a one-crossing solution 
while figure (ii) does the same for a two-crossing solution. The third figure illustrates the 

solution for each p. The solution of p is ^^^^j^sinh \J~^^'r] in fixed AdS space. Thus, a graph 

of hyperbolic sine type function near the vacuum states indicates AdS space. The maximum 
number Umax is 4 for these parameters Uo and k. 

Figure [8] shows the collective variation of terms in equations of motion, phase diagrams, 
and the diagram of energy density between AdS-AdS degenerate vacuum states. Figure (i) 
corresponds to a one-crossing solution while figure (ii) corresponds to a two-crossing solution. 
In the third and sixth figures, the dented region between the peaks of the Euclidean energy is 
due to the integration of variables in AdS space. The volume energy density and the Euclidean 
energy are both negative in AdS region and become positive in the valley of the inverted 
potential, i.e. the dS region. The initial and the final region of p' in the first and fourth figures 
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Figure 5: The numerical solutions between flat-flat degenerate vacua. 

exhibit a hyperbolic cosine type function, because the solution near the vacuum states indicates 
AdS space. In Figs. [7] and [HI two solutions with n = 1 and n = 2, are illustrated and all other 
solutions are omitted since the general behaviors for these cases are similar to those in dS space. 

The second term of the first equation in Eq. (jlj) contributes significantly to obtain solutions. 
In other words, the sign of the term needs to be properly changed to keep harmony with the 
sign of the other terms. This adequate change of sign is possible only if the gravity is taken 
into account. As a result of tunneling, geometries of solutions become finite with Z2 symmetry. 
Each space surrounded by the wall can be either dS, flat, or AdS space. Thus the geometries 
of solutions in flat and AdS space become quite different from those before tunneling. 

We now examine the tunneling probability for the solutions. The tunneling rate can be 
evaluated in terms of the classical configuration and represented as e"^*^ in this approximation, 
where the leading semiclassical exponent AS = S'^^ — S^^ is the difference between the Euclidean 
action corresponding to a classical solution S'^^ and the background action S''^. When one 
considers the usual bounce solution, the outside geometry of the bounce solution does not change 
ever after tunneling. Hence the contribution from outside is canceled. The only contribution 
to AS* comes from the part inside and the wall leading to the finiteness of AS. 

On the other hand, if one considers a single instanton solution with a finite size, there exists 
a subtle problem. There is no such subtlety for the dS case. For this case, both the actions 
for solutions with a finite size of the lens type geometry and background with finite size S'^ 
are negative, i.e. S"^* < and S''^ < 0, while the net difference is positive AS' > giving 
rise to the probability of less than 1. However, there is a problem in defining the probability 
for the flat and AdS space. For flat space the Euclidean action for solution is a negative one, 
i.e. 5''^'* < while the background S^^ = 0. Therefore the net difference is negative AS* < 
making the naive probability greater than 1. The negative value of the exponent is due to 
the contribution for the valley of the inverted potential, i.e. the dS region, which is related to 
the fact that the Euclidean action for Einstein gravity is not bounded from below [25]. For 
AdS space the Euclidean action of solutions is a negative one S'^''^ < 0, whereas the background 
5**^ = +00. The net difference has got the negative divergence AS" = —00. This divergence 
is due to the contribution coming from the infinite space with nonzero vacuum energy. It is 
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Figure 6: Variation of terms in equations of motion, phase diagrams, and the diagram for energy density 
between flat-flat degenerate vacua. 



not clear how to interpret this point in physical terms in the present work. We can make AS* 
finite by introducing a cutoff. The effect of the cutoff, as the size increases, was discussed in 
Ref. [7]. In this work, we will concentrate only on the relative probability. We then only need 
to compare the actions for oscillating solutions with the action of a one-crossing instanton, i.e. 

foes Cbg\ ( Qcs Cbg\ ocs ocs 

Wn ) Wl ) — '-'n '-'l ■ 

Table [1] shows the actions for all of the solutions. We numerically obtained the actions for 
oscillating solutions from the data. In the integration, we employ an additional normalization 
process as we use dimensionless quantities for our numerical calculation. Thus we have nor- 
malized ASn to be the difference between the action of the n-crossing solution and that of a 



one-crossing solution divided by a 



Si 



as follows: 



A5„ 



Sn — Si 



Si 



Si 



(16) 



where Sn denotes the action for the ra-crossing solution. Our results are depicted in Table [H 

in which (^k,Uo^ are taken to be (0.04,0.5), (0.2,0), and (0.4,-0.02) respectively. Our results 

show that, as predicted, the transition amplitude is suppressed with an increase in oscillation 
number. 



13 



20 40 10 20 "o 15 30 45 

r| 11 11 



Figure 7: The numerical solutions between AdS-AdS degenerate vacua. 
Table 1: A5„ for Oscillating solutions between dS-dS, flat-flat, and AdS-AdS degenerate vacua. 



Oscillation number 


dS-dS case(A^„)(e-^^") 


flat-fiat (A5„)(e-^^") 


AdS-AdS(A5„,)(e-^^") 


1 


-544000(0)(1) 


-403752(0)(1) 


-152000(0)(1) 


2 


-511000(0.061)(0.941) 


-161000(0.601)(0.548) 


-39300(0.903) (0.405) 


3 


-490000(0.099)(0.906) 


-114726(0. 716)(0.489) 


-30234(0.925) (0.396) 


4 


-479000(0. 119)(0.887) 


-101000(0. 750)(0.472) 


-28384(0.930) (0.395) 


5 


-475000(0. 127)(0.881) 


-96000(0.762)(0.467) 




6 


-473000(0. 131)(0.878) 


-94700(0.765) (0.465) 





4 Properties of oscillating instanton solutions between 
the degenerate vacua 

In the previous section, we obtained various properties of the oscillating solutions between dS- 
dS, flat-flat, and AdS-AdS degenerate vacua. We make an attempt to obtain the phase space of 
these solutions in terms of the parameters k and klJo- The parameter k is the ratio between the 
gravitational constant or Planck mass and the mass scale in the theory, k = = ^ly, whereas 

the parameter kUo is related to the cosmological constant A//i^. Many questions naturally arise 
at this stage. The first one is how many of the oscillating solutions are allowed for given k and 
Uo- The second one is whether the number of oscillations depends on these parameters or not. 
What is the whole phase space of solutions according to the parameter regions, if the number 
depends on the parameters? The third one deals with the different behaviors of solutions among 
dS-dS, flat-flat, and AdS-AdS cases. To answer these questions, we will try to figure out the 
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phase space of solutions in terms of R and kUq- For k <^ 1, the effect of gravity is neghgible. 
However, when k approaches to order one value, the effect of gravity becomes important. To 
obtain the phase space of solutions, we collect data of f/o, which makes a difference in the 
maximum number of oscillations, Umaxi for a given k. We then obtain a specific Uo which 
determines the minimum number of oscillations, rimm- Finally, we employ the method of least 
squares [2S] to obtain the relationship between the two parameters in the given data sets of k 
and kUo- Oscillating numbers appear to be linearly related to the parameters. 

Figure in] shows the behavior of oscillating solutions in terms of k and klJo- The value of k is 
limited in the range ^ k ^ 1. At the point k = 0, or Mp = oo, gravity is switched off. The y 
axis represents no gravity. We expect that there exists no solution with 0(4) symmetry when 
gravity is turned off. In each of the figure parts, we use the notation {rimm, nmax), where rimm 
represents the minimum number of oscillations and rimax f he maximum number of oscillations 
in the given parameter range. For example, {rimini nmax) = (l? 1) nieans that the minimum 
number of oscillations is 1 and the maximum number of oscillations is 1. Figures (i) and (ii) 
illustrate the number of oscillations for the case of dS-dS degenerate vacua. In fig. (i), there is 
a zone representing no solution in the upper right region. When the parameters k or Uo are 
increased, the evolution period of rjmax is diminished in general, p approaches pmax before the 
field arrives at the other vacuum state in the inverted potential. Thus, there is no solution 
because of the short evolution period. In other words, the instanton solution can not fit inside 
the Euclideanized dS background in the strong gravity limit. When the parameters k or Uo 
are decreased, both the evolution period of r/max and the maximum number of oscillations get 
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Figure 9: Oscillating behavior depending on k and kUo- Figure(i) shows the oscillating properties of dS-dS 
degenerate case and figure(ii) shows small values of parameters once more. Figure(iii) is a AdS-AdS degenerate 
case, and figure(iv) shows the total parameter space, we have searched. 



increased. Figure (ii) illustrates the zone with small value of parameters in fig. (i). The upper 
zone allows rimin = 1? whereas the lower zone allows rimm = 2 divided by the solid line. Figure 
(iii) illustrates the behavior of the oscillating solutions in the case of AdS-AdS degenerate vacua. 
In this case, {nmin-,nmax) = (1, 1) is not allowed. If there is a tunneling solution, then there are 
oscillating solutions found together. This is a different property as compared to the dS-dS case. 
In addition to this, there is a no-solution parameter region corresponding to regions for which 
Lfo < —0.125 in the case of AdS-AdS degenerate vacua as there is no dS region of potential [7]. 
In this case, the change of n^m is more prominent. Even at k = 1, with strong gravity, the 
change of rimin occurs. This is another difference between the dS-dS and the AdS-AdS cases. 
Whole parameter space of R and kUo is shown in the fig. (iv). In this figure, we can see that 
the changes of rimin and Umax in the dS-dS case are continuously connected to those in the 
AdS-AdS case. The x axis, Uo = 0, represents the fiat-fiat case. In the fiat-fiat case, there 
exists a solution for all of k except for k = 0. 

Figure [TU] shows the schematic diagram for the phase space of all solutions including yet 
another type of solution and the number of oscillating solutions with different k values. The 
left figure has k = line indicating no gravity effect. There is a zone representing no solution in 
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Figure 10: Illustrations of the phase space of solutions and oscillating properties. Left figure shows a whole 
solution space, and right figure shows number of oscillating solutions changing on given k. 



the upper right region in the case of dS. The dS region has positive Uo and the AdS region has 
negative Uo divided by the hne, Uo = 0, representing the flat case. In the middle area including 
the flat case, rimin and and kUo decrease. The tendency is indicated as the 

painted arrows. The inchned line in the AdS region represents Uo = —0.125 or Utop = 0. In the 
lower left region, there exists another type of solution. The figure on the right shows how n^m 
and Umax are changed in terms of kUo and k. As we can see from the figure, rimax and n^m 
increase as kUo decreases. 

Figure [11] shows another type of solutions in Fig. [TOl The solutions represent tunneling 
from the top of the potential, a point of an unstable equilibrium, to the local vacuum state 
instead of rolling down the potential. The solutions can be of the same kind of a bubble solution 
describing tunneling without a barrier [311132] . In figures (i), we take k = 0.5 and Uo = —0.125. 
The inside geometry is AdS and that of the outside is flat. In figures (ii), we take k = 0.5 and 
Uo = —1. The inside geometry is AdS and that of the outside is AdS. 



5 Summary and Discussions 

We have studied oscillating instanton solutions of a self-gravitating scalar field between de- 
generate vacua. We obtained numerical solutions in dS background. Basically, our method 
for obtaining the domain wall or brane world- like object [331134] is based on the instanton- 
induced theory rather than the kink-induced theory. Our approach is related to the question: 
How we can make spacetime including the domain wall in dS and flat space? If the thin-wall 
approximation scheme is allowed in our work and the mechanism is applied to solutions in 
higher-dimensional theories, our one-crossing solutions can be interpreted as the mechanism 
providing nucleation of the domain wall or braneworld in instanton-induced theory. Because 
our oscillating solutions have a thick wall with varying energy density, our oscillating solutions 
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Figure 11: The solutions describing tunneling without a barrier. The upper two figures show the solution of 
$ and the lower two figures show the solution of p. 

can be interpreted as the mechanism providing nucleation of the thick wall for topological 
inflation [35|[36]. We add that Z2 invariant solutions also exist in flat or AdS background, 
though the physical signiflcance is not clear. Furthermore, we constructed the phase space of 
all our solutions including the number of oscillations. As a by-product, we obtained the solution 
describing tunneling without a barrier. 

In Sec. II, we analyzed the boundary conditions of our problem. There are two kinds of 
conditions. We imposed the boundary conditions as the initial value problem and imposed 
additional conditions implicitly. The initial value $0 is obtained by employing the undershoot- 
overshoot procedure. To avoid a singular solution at f/ — )■ fjmax in Eq. flTU]) and demand the Z2 
symmetry, the conditions d^/drj — )■ and p as fj ^ fjmax were needed. 

In Sec. Ill, we obtained the numerical solution of oscillating instantons. In particular, we 
performed the numerical work with more detail for the cases in dS background. The number 
of oscillations increases as the initial point, ^{fjinuiai) = ^o, moves away from the vacuum 
state. We can see that the size of the geometry of solutions decreases as the number of crossing 
increases because the period of the evolution parameter fj^ax decreases as the starting point 
moves away from the vacuum state. As expected, this type of solutions is possible only if gravity 
is taken into account. The maximum number of oscillations is determined by the parameters k 
and Uo observed in Ref. [23]. The solutions are of two types. One type is for situations in which 
different parts of spacetime are in different vacua. The solutions representing tunneling which 
go from the left vacuum state to the right vacuum state have a nonzero topological charge. The 
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other type is for situations in which different parts of spacetime are in the same vacuum state. 
Others representing the solutions which go back to the starting point after oscillations have 
zero topological charge. We can see that the sign change of p' from positive to negative occurs 
at the half period of fj due to the Z2 symmetry. We have studied the behavior of the solutions 
in the -$'(^) plane using the phase diagram method. Figures [3]|^i), (iii), and (v) have 
symmetry about the y axis, while figures [3](ii), (iv), and (vi) have symmetry about the x axis. 
The maximum value of $' decreases as the number of crossing increases. We have checked the 
energy density diagrams for each solutions. The peaks of the volume energy density broaden 
in range near Utop as the number of crossing increases. The shape of peaks of the Euclidean 
energy becomes smooth and broadens as the number of crossing increases. 

We have examined the tunneling probabilities for the solutions. The tunneling rate can 
be evaluated in terms of the classical configuration and can be represented as e"^*^ in this 
approximation, where the leading semiclassical exponent AS* is the difference between the 
Euclidean action corresponding to a classical solution and the background action. When one 
considers the usual bounce solution, the outside geometry of the bounce solution does not 
change even after tunneling. Hence the contribution from the outside gets canceled. The only 
contribution to AS* stems from the inside part and the wall leading to the finiteness of AS*. 

On the other hand, if one considers a single instanton solution with finite size, there exist 
a subtle problem. There is no such subtlety for the dS case. However, there is a problem in 
defining the probability for the fiat and the AdS background. The Euclidean action for the fiat 
space case is a negative one, AS < making the naive probability greater than 1. The Euclidean 
action for the case of AdS space AS* — ?■ — 00. This negative divergence of the exponent is due 
to the contribution coming from the infinite space with nonzero vacuum energy. It is yet not 
clear how to interpret this point correctly in the present work. Of course, we can make AS 
finite by introducing a cutoff. The effect of the cutoff, as the size increases, was discussed in 
Ref. [7]. In the present work, we concentrate mainly on the relative probability. We only need 
to compare the action of oscillating solutions with that of a one-crossing instanton. Our results 
show that, as predicted, the transition amplitude is suppressed with an increase in the number 
of oscillations. One may concern the contribution from the boundary term [2H]- We expect 
that the point fj = fjmax in the opposite side after the tunneling process is smooth due to Z2 
symmetry and therefore the boundary term does not contribute to the action [7]. 

In Sec. IV, we obtained the phase space of solutions in terms of the parameters k and kUo- 
To make the phase space of solutions, we collected data on Uo, which makes the difference 
of the maximum number of oscillations, rimax, for a given k. We then obtained a specific Uo 
which determines the minimum number of oscillations, rimin- Finally, we employed the method 
of least squares [2B] to obtain the relationship between two parameters in given data sets of 
k and kUo- Oscillating numbers appeared to be linearly related to the parameters. Figure M 
shows the behaviors of oscillating solutions in terms of k and klfo- The value of k is limited 
to ^ K ^ 1. At the point k = 0, or Mp = 00, gravity is switched off. The y axis represents 
no gravity. We expect that there is no solution with 0(4) symmetry when gravity is switched 
off. When the parameter k or Uq increased, the evolution period of "fj^ax diminished in general. 
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When the parameter k or Uo decreased, the evolution period of fjmax and the maximum number 
of oscillations increased. There is a no-solution parameter region corresponding to regions in 
which Uo < —0.125 in the case of AdS-AdS degenerate vacua because there is no dS region of 
potential [7]. We can see the changes oirimin and rimax in the dS-dS are continuously connected 
to those in the AdS-AdS. The x axis, f/o = 0, represents the fiat-fiat case. In the flat-flat case, 
there exists a solution for all of k except when k = 0. Figure [10] shows the schematic diagram 
of the phase space of all solutions including another type solution and the number of oscillating 
solutions with different ks. 

As a result of this tunneling, a finite-sized geometry with Z2 symmetry is obtained. Our 
mechanism for making the domain wall or braneworld-like object is different from the ordinary 
formation mechanism of the domain wall because our solutions are instanton solutions rather 
than soliton solutions. In other words, our solutions can be interpreted as solutions describing 
an instanton-induced domain wall rather than a kink-induced domain wall or braneworld-like 
object. Domain walls can form in any model having a spontaneously broken discrete symmetry. 
An inertial observer sees the domain wall accelerating away with a specific acceleration. The 
domain wall has repulsive gravitational fields [nZllSB]- The thickness of the domain wall in 
flat spacetime can be estimated by a balance between the potential energy and the gradient 
energy. When the thickness of the domain wall is greater than or equal to the horizon size 
corresponding to the vacuum energy in the interior of the domain wall, topological inflation 
can occur. The scalar field stays near the top of the potential at the core. This potential 
energy serves as a vacuum energy in a similar way to the slow-rollover inflationary models. 
This topological inflation does not require fine-tuning of the initial conditions and is eternal 
even at the classical level due to the topological reason. Our oscillating instanton solutions 
can be interpreted as mechanism providing the nucleation of the thick wall for the topological 
inflation. The wrinkles representing the variation of the volume energy density in the wall may 
be interpreted as density perturbations in the inflating region. In this work, inflating regions 
described by the oscillating solutions and density perturbations described by the variation of 
energy density can occur simultaneously. Furthermore, oscillating bounce solutions also have 
the thick wall. Thus we expect that (non-)topological inflation can be made by oscillating 
bounce solutions. 

To obtain the dynamics of the solutions, we should apply the analytic continuation from 
Euclidean to Lorentzian signature. If the wall is thin, we can employ the Israel junction 
condition [TJIMllin] or the method [lOlIlT]. If the wall is thick, the double-null simulation may 
be more relevant to the dynamics of the solutions |42] . 

A direction for future research could be whether or not our solutions obtained using the 
instanton-induced theory can be extended to theories with gauge fields or various other di- 
mensions and whether or not the topological charge of both instanton and bounce solutions, 
including oscillating one, can be well-defined and conserved for a self-gravitating scalar field in 
various dimensions. It would also be interesting to examine if this toy model with the variation 
of energy density can provide a proper inflationary scenario. In Ref. [13], the creation of a 
charged black hole pair separated by a thin domain wall and the dynamics of the domain wall 
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were studied in the cosmological context. It would be interesting to study the properties and 
evolution of topological inflation with a magnetic fleld. 
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